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ABSTRACT 

Sterile neutrinos may be one of the best Warm Dark Matter candidates we have today. 
Both lower and upper bounds on the mass of the sterile neutrino come from astro- 
nomical observations. We show that the proper inclusion of the neutrino momentum 
distribution reduces the allowed region to be 2.6 keV < m < 5 keV for the simplest 
models. A search for a spectral line with E = to/2 is thus more interesting than ever 
before. 
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Astrophysics provides an increasing amount of independent 
indications that the dark matter of the universe is warm, 
so that the small-scale fluctuations are damped out by free 
streaming. This is most easily achieved by giving a keV mass 
to the DM particle, in which case the preferred candidate is 
the sterile neutrino. Support for warm dark matter (WDM) 
comes from simulations of the number of satellite galaxies 
(Colin, Avila-Reese & Valenzuela 2000) and of disk galaxy 
formation without the need for stellar feedback (Sommer- 
Larsen & Dolgov 2000), which both find that a DM particle 
mass of about 1 keV is optimal: a significantly larger mass 
has little impact on galaxy formation, and a significantly 
smaller mass would lead to the well known difficulties faced 
by hot dark matter. A quantitative lower limit on the can- 
didate WDM particle mass is inferred from the existence 
of a massive black hole at large redshift (Barkana, Haiman 
& Ostriker 2001) and the requirement of sufficiently early 
galaxy formation to account for reionization of the universe 
and the observed Ly-a forest properties (Narayanan et al. 
2000), constraining the DM mass to be larger than 0.75 keV. 
A recent discussion of x-ray emission from decays of sterile 
neutrinos (Abazajian, Fuller & Tucker 2001) has imposed an 
upper limit of about 5 keV on the neutrino mass. Here we 
discuss a reinterpretation of these bounds on neutrino mass, 
and demonstrate that the proper inclusion of the neutrino 
momentum, arising from the specific production tempera- 
ture, reduces the allowed sterile neutrino WDM mass to be 
in the range 2.6 keV < m < 5 keV. 



2 WDM PARTICLE DECOUPLING 

All of these studies (Colin et al. 2000; Sommer-Larsen & Dol- 
gov 2000; Barkana et al. 2001; Narayanan et al. 2000) are 
based on the mass-dependent cut-off on small scales, pro- 
duced by free-streaming. In the previously cited studies, a 
"conventional" WDM model was considered for the under- 
lying particle physics (see Bode, Ostriker & Turok (2001) 
and Sommer-Larsen & Dolgov (2000) for recent overviews 
of such particle models). In such conventional WDM mod- 
els, the particles decouple in the early Universe at higher 
temperatures than do massless neutrinos. Therefore they do 
not share the entropy release from the successive particle an- 
nihilations. Since they were relativistic at decoupling, their 
distribution function in momentum space is subsequently 
that of a massless fermion, but with a temperature, Tw, 
which is given today by 



94eV\ 1/3 



inw 



(1) 



where T„ « 1.946 K, H = 100 /ikms _1 Mpc _1 , tt w is the 
present energy density of WDM in units of the critical den- 
sity, and mw is the WDM mass. The needed entropy release 
in these conventional models is much bigger than allowed 
in the standard model, and such WDM candidates should 
therefore have decoupled before a larger gauge group breaks 
down. 

Now a very natural candidate for the WDM particle is 
a massive sterile neutrino mixed with an ordinary neutrino 
(Dodelson & Widrow 1994; Colombi, Dodelson & Widrow 
1996; Shi & Fuller 1999; Dolgov & Hansen 2001; Abazajian, 
Fuller & Patel 2001). Since the mixing angle is temperature 
dependent (Notzold & Raffelt 1988) (and for small vacuum 
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mixing angle, sin 2 26 ~ 10~ 7 ), only a small amount of these 
heavy neutrinos, relative to ordinary active neutrinos, can 
be produced at high temperatures. The distribution function 
of sterile neutrinos is, to a fair approximation, characterized 
by the temperature of massless neutrinos, but smaller by a 
factor X. For a specific choice of mw and Slw today the 
value of X can be found from 



x = n, 



,2 /94eV\ 

rh ~ 10 

V mw ) 



(2) 

mw 

therefore the two models produce the same contribution to 
ptot of WDM particles today if 

V i —2- I . (3) 



However, WDM particles with mw have a different distribu- 
tion function in these two models, and their free-streaming 
effect is not equivalent. The effect on large scale structure 
was first discussed in detail by Colombi et al. (1996). We 
will use conventional WDM (cWDM) when referring to the 
first case, and sterile neutrino WDM (sWDM) for the second 
one. 



3 LOWER BOUNDS 

The two neutrino models are easily included in a Boltzman 
code, in order to compute the present matter power spec- 
trum P(k). Using the code cmbfast (Seljak & Zaldarriaga 
1996), we have found analytical fits for the transfer func- 
tions, T(k), relating the power spectrum in the WDM to 
the CDM scenario 



T 2 (k) 



P w (k) 
P CDM {k) 



for Ot, 



ficz 



(4) 



where P w is the power spectrum for cWDM, and a similar 
expression with P h ' (k) for the sWDM model. These transfer 
functions, which essentially reflect the free streaming cut-off, 
have the form 

T(k)= [l + {akrY 5/ \ (5) 

where k is the wavenumber in units /iMpc -1 , v — 1.12, and 
a depends on the cosmological parameters as 



a ( nw \ b ( h V ( mw V 

V0.3/ Vo.65/ UoOeV/ 



(6) 



Numerically we find for cWDM, A = 1.07, 6 = 0.11, c = 1.20 
and d = —1.11 in good agreement with Bode et al. (2001). 
For the sWDM one can derive similar numbers by noting 
that the mass in the cWDM case differs by (Tw /T VQ ). 
This means that if we have a dependence Qwh c m w for the 
sWDM case in eq. and a dependence 0. w h c m w for the 
cWDM case, then one finds 



Q, w h c m w = 



b'-d'/3 ,c'-2d'/3 d'+d'/3 



(7) 



which is solved by b' = b + d/4, d' = 3d/4 and c' = c + d/2, 
in good agreement with what we found numerically, by ex- 
plicitly changing the massive neutrino phase-space distri- 
bution function as done by Lesgourgues & Pastor (1999). 
To be very explicit, this means that for a given cut-off 
scale of the power spectrum one can find the correspond- 
ing mass of the elementary particle, and the mass will differ 



in the two cases. E.g. if for cWDM one finds mw = 0.75 
keV, then this corresponds in the sWDM case to mw = 
(fiw-fr 2 94eV/0.75keV) 1/3 ■ 0.75keV « 2.6 keV, when us- 
ing h = 0.7 and Qw = 0.4. In other words, if one believes 
that sterile neutrinos indeed constitute the dark matter, 
and they are produced as described in (Dodelson & Widrow 
1994; Colombi, Dodelson & Widrow 1996; Dolgov & Hansen 
2001; Abazajian, Fuller & Patel 2001), then the bounds ob- 
tained by Barkana et al. (2001) and Narayanan et al. (2000), 
mw > 0.75 keV, should really be multiplied by a factor 3.4, 
and the lower bound on sterile neutrinos as dark matter is 
thus about 2.6 keV. 

This lower bound may be subject to minor correc- 
tions. First, the temperature of the sterile neutrinos is really 
slightly lower than the active neutrino temperature, since 
the sterile neutrinos are being produced while the muons 
are still present in the Universe, T ps 130 MeV (Langacker 
1989; Kainulainen 1990; Barbieri & Dolgov 1990, 1991). Sim- 
ilarly, for the neutrino states being produced above the QCD 
phase transition, one must also take into account the quark 
degrees of freedom (Abazajian et al. 2001a). Another effect 
arises from the fact that the produced sterile neutrino spec- 
trum is not exactly thermal, but slightly warmer in the sense 
that the higher momentum part is more populated than the 
lower momentum part (Dolgov & Hansen 2001; Abazajian 
et al. 2001b). Furthermore, the factor of 3-4 found above de- 
pends on the specific values of flw and h, and can therefore 
change slightly. It is also worth noting, that if the sterile neu- 
trinos are produced resonantly (Shi & Fuller 1999) , through 
a pre-existing lepton asymmetry, then the upper limit on 5 
keV may weaken substantially (Abazajian et al. 2001b). 



4 UPPER BOUNDS AND DISCUSSION 

Sterile neutrinos in the keV mass range have a decay time 
that is of cosmological interest. Recently a very interesting 
paper appeared (Abazajian et al. 2001b), where the signa- 
ture from decaying sterile neutrinos in galaxies and clusters 
of galaxies was studied in detail 0. A bound m < 5 keV was 
derived, using the relation between mass and mixing angle 
obtained in (Abazajian et al. 2001a). We thus conclude that 
a sterile neutrino as WDM must lie in the mass range 2.6 
keV < m < 5 keV, and it is therefore more interesting than 
ever before to search for a spectral line with energy E = m/2 
from the decay v s — > v a + 7. 

If future searches for a spectral line from the sterile neu- 
trino decay should give a negative result, then one must find 
new WDM candidates. An interesting possibility is an ac- 
tive neutrino, which may never reach thermal equilibrium if 
the reheat temperature at the end of inflation is low enough 
(Giudice et al. 2001). Such a scenario demands a reheat tem- 
perature of a few MeV, and the resulting neutrino distri- 



* One could imagine a slight change of strategy in the analy- 
sis of (Abazajian et al. 2001b), namely to consider regions with 
large concentration of dark matter, but with little baryonic mat- 
ter. Such "dark blobs" may have been observed by inverting the 
matter distribution in clusters of galaxies from weak lensing, and 
comparing with the baryonic matter inferred from optical obser- 
vations of the cluster (Clowe et al. 2000), but the significance of 
such blobs is far from being established. 
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bution function is also warmer than a conventional WDM 
candidate. Such a solution does, however, not appear too 
natural in view of the recent neutrino data (Ahmad et al. 
2001; Toshito et al. 2001), indicating that all the active neu- 
trino masses are sub eV. Several WDM candidates have dif- 
ferent bounds from what is presented here, e.g. gravitinos 
(Pagels & Primack 1982; Borgani, Masiero & Yamaguchi 
1996; Kawasaki, Sugiyama & Yanagida 1997) produced in 
the very early Universe; or sterile neutrinos if there should 
exist an initial lepton asymmetry (Shi & Fuller 1999). Even 
more traditional DM candidates (like a neutralino) can dis- 
guise themselves as WDM, namely if they have scattering 
cross section with photons or neutrinos. Specifically, for a 
100 GeV DM particle with scattering cross section about 
(tdm- 7 ~ 10 _30 cm 2 one obtains a reduction of the mat- 
ter power spectrum corresponding to a conventional WDM 
candidate with mass about 1 keV (Bcehm et al. 2001). 

In conclusion, we have shown that the proper inclusion 
of the neutrino momentum distribution changes the lower 
bound allowed for the simplest models of sterile neutrinos 
as WDM, and the resulting allowed region thus becomes 2.6 
keV < m < 5 keV. 
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